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Abstract 

A generic PT— symmetric Hamiltonian is assumed tridiagonalized and truncated to 
N < oo dimensions, H jjicham model) ^ g^jj -^.g yp_(;Jown symmetrized special 
cases with J = [N/2] real couplings are considered, H^^^'^'^^ model) _^ H^^\ Using 
symbolic manipulation and extrapolation techniques we find out that in the strongly 
non-Hermitian regime the secular equation gets partially factorized at all A^. This 
enables us to reveal a fine-tuned alignment of the dominant couplings implying an 
asymptotically sharply spiked shape of the boundary of the J— dimensional quasi- 
Hermiticity domain D^^) in which all the spectrum of energies E^^^^ remains real and 
observable. 
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1 Introduction 



Unexpectedly often, many general theoretical considerations as well as practical phe- 
nomenological applications of quantum mechanics rely on the exceptionally friendly 
mathematics connected with the one-dimensional harmonic oscillator H^^^\ In par- 
ticular, the equidistance of its energies (say, = 1, 3, 5, ... in suitable units) proves 
fairly favorable in perturbation theory where it simplifies the practical calculations 
of the spectra of the various anharmonic oscillators 

■j^(AHO) _ jj{HO) _j_ (perturbation) ^j^^ 

Even when one restricts attention to the mere finite- dimensional perturbations writ- 
ten in terms of the eigenvectors | m^^^^) of H^^^^ itself, 

N 

m, n — 1 

the finite-dimensional matrix example ([I]) + ([2]) may produce a lot of theoretical 
inspiration as sampled, e.g., in Chapter Two of the Kato's classical book on per- 
turbation theory where many low- dimensional Hermitian as well as non-Hermitian 
sample matrices were considered 

The parallel, purely phenomenological inspiration by the truncated eq. ([2]) need 
not be less exciting. For example, the authors of an explicit numerical exercise 
[2] felt inspired by the recent growth of popularity of the PT— symmetric f5] and 
pseudo-Hermitian P] models and analyzed their strongly-perturbed sample with 

W = i (^(HO) I ^3 ^iHO) ] I ^(HO)^ ^ m, n = 1, 2, . . . , iV. (3) 

They found out that as long as H^^'^^ = x^+p^ is positive, a very good approximation 
of the nonlinear-differential-equation N = oo results is achieved via the series of 
truncated, purely algebraic and linear N < oo eqs. ([2]). 

In our recent papers [5] and [6] we made one more step. Assuming that one 
could annihilate any far-off-diagonal element of any given matrix H^^^^^ via a finite 
sequence of elementary Jacobi rotations [7] we restricted our attention to the mere 
"irreducible" , tridiagonal anharmonic-oscillator-type models 

N 

^icha^n) ^ ^iHO) ^ ^ | m^^^^) W^ifj (n^^^) | . (4) 

m, 71 — 1 
\m — n\ — 1 
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In the PT— symmetric scenario (guaranteed, in the normahzation accepted in [5] , by 
the mere antisymmetry W, 



(N) 

m,m+l 



-W 



(N) 



m+l,m 



Qm among the real matrix elements) 
we shifted the origin of the energy scale and arrived at the A^— dimensional and 



tridiagonal "chain-model" matrices {m^^'^^ \ if 



n 



jjichain) 



-N 

-9i 
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(HO) 






\ — TT{chain) 
I m,n ' 



9N-2 

N -3 



9n~i 



-g^.i N-1 



(5) 



For the sake of simplicity we decided to pay attention solely to the up-down-symmetric 
special cases H^'^'> of if ('='*<^™) where we choose 



gN-k = 9k>^ 



A; = 1,2,..., J 



(6) 



at both the even N = 2J and the odd = 2 J + 1. 

In the subsection 12. II of section [2] below we shall start our present continuation of 
the latter study of the toy model ([5]) + ([6]) by a brief review of the results of paper [5] . 
We point out there that the model H^^^ itself has been introduced as partially 
tractable by an algebraic, symbolic-manipulation-based non-numerical extrapolation 
method. We remind the readers that at any dimension = 2J or = 2J + 1, 
the spectrum {E^^} remains real and observable inside a J— dimensional domain 
T) = pC^) of the matrix elements which is compact and all contained inside a bigger 
domain S'-^'^ defined by the following elementary inequality [S], 



-N 




N 
N 



2 J, 

2J + 1. 



(7) 



There exists just a finite set of the "maximal-coupling" intersections of the two sur- 
faces, i.e., of the (J — 1)— dimensional boundaries &D^^^ and dS^^\ We succeeded in 
determining the coordinates of these points (called, in ref. [5], extremely exceptional 
points, EEP) in closed form. 



niEEP) 



n] 



n 



1,2, 



J. 



In the subsection 12.21 we shall extend the review by adding some empirical observa- 
tions published in our most recent numerical study ^ and concerning the behaviour 
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of the energies predominantly, far off the EEP extremes. In particular, we 

shall recollect there the lucky guess of the ansatz 



9n = ^{l-Ut)) , Ut)=t + t' + ... + t'-' + G^t' (9) 



which extrapolates, to all J, the rigorous fine-tuning rule derived, in ref. [5], at J = 2. 
In the numerical context of ref. |6j it merely served as a bookkeeping tool in our 
experiments with the various choices of the rescaled couplings Gn- In what follows 
we intend to describe several much deeper and more far-reaching consequences of 
this type of an ansatz. 

In the preliminary steps made in section[3]we shall start the analysis of our bound- 
state problem by the direct, brute- force algebraic solution of its secular equations 



at = 4 (subsection 13.11) . at = 5 (subsection 13.21) and at = 6 (subsection 13.31) . 
In the subsequent section H] we shall change the variables in the manner prescribed 
by eq. iQ. This will lead to the much more compact strong-coupling leading-order 
formulae at the two sample dimensions = 4 (subsection l4.1l) and A^ = 5 (subsection 
14.21) . Certain indications of the possibility of a successful extrapolation of these 
formulae with respect to the dimension A^ will follow in section [3] sampling A^ = 6 
(subsection 15.11) and N = 7 (subsection 15. 2p . 

The climax of our present paper comes in sections [6] (where we extend the above 
results to all the even dimensions A^ = 2 J) and[7| (where the parallel extrapolations 
are outlined in the case of any odd A^ = 2J + 1). The remaining text are just 
discussions (section [8]) and summary (section [9]) which emphasize that our present, 
perturbation-theory-simulating attention paid to the interval of small t fills in fact 
the gap between the algebraic t = approach of paper and its numerical large— t 
pendant of letter [B] . 

2 A brief review of the state of art 

Our continuing interest in the family ([5]) has several reasons. Firstly, in possible 
applications, special cases of H^^^ could play the role of a non-standard spin model 
(cf. their N = 2 samples in refs. P, [H]) or of a Hamiltonian of a system where 
a "de-frozen", new degree of freedom can emerge (with A^ = 3 see [lOj). At any 




(10) 



4 



higher dimension N, every is, using the language of the review paper [TT] . 

quasi-Hermitian and, therefore, ehgible, say, as a Hamiltonian of a quantum chain. 
Inside V^^^ and in a suitably specified physical Hilbert space H^^\ all our models 
H^^^ obey all the postulates of Quantum Mechanics in a way illustrated, at = 2, 
in [9]. Due to their finite matrix form, all of them appear particularly suitable 
for an illustration of some subtleties of perturbation theory, especially in the weak- 
coupling dynamical regime where all their off-diagonal elements Qn remain small [1] . 
In the sense explained in [5], all of our toy Hamiltonians H^^^ are also interesting 
as PT— symmetric [3l [12] and/or parity-pseudo-Hermitian [3] candidates for opera- 
tors of observables which are capable of exhibiting multiple confluences of their (in 
the Kato's language) "real exceptional points" [HI |13l IHj. Last but not least, our 
non-Hermitian matrices H^^^ are, through the above-mentioned variational and per- 
turbative considerations [15j, directly and closely related to the popular differential 
PT— symmetric quantum Hamiltonians of Bender et al PJ] et al [TB] - |21j . 



2.1 Observations made in paper |[5| 

In "paper I" [5] we felt inspired by the exceptional transparency of the geometry of the 
"physical" domains of quasi-Hermiticity T'(^) in the strong-coupling regime and at 
the smallest dimensions A^. In such a setting, the first non-numerical result of ref. [5] 
was that all the eigenenergies become complex whenever the anharmonicity becomes 
sufficiently strong. Such a type of observation was interpreted as very important 
because the crucial point of making any PT— symmetric Hamiltonian H "physical" 
(i.e., responsible for observable and stable bound states) lies in the specification of 
the "allowed" range V = T>{H) of its free parameters. 

In paper I we showed that the non-nnmenca/ construction of 'D{H) proves feasible 
in an EEP, "maximal-coupling" limit H = Hj^^l,p^ of our tridiagonal matrices H^^^ 
at all their dimensions N. The proof was based on an elementary observation that for 
every individual model H^^^ with a fixed dimension N, the spectrum is determined 
by the polynomial secular equation for the squared energy s = E'^, 

det [i/ - E /(^)] =s'- Pf ) {g,,..., gj) s^'^ + P^^^ {g„ . . . , g.j) s'-'- 

-... + (-l)^pf)(^i,...,^,)=0. (11) 
We recollected that due to the polynomiality of the latter equation, the sum of the 
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physical (i.e., nonnegative) roots Sj must be equal to the first coefficient, 

Pi''\gi, ...,gj) = si + S2 + ... + sj>0. (12) 

As a consequence, the closure of the domain p(^) must lie inside the closure of an- 
other, bigger domain S^'^^ which is defined, much more simply, by the upper estimate 
ffT^ . The shape of the surface of S^^^ is a hyperelhpsoid or a hypersphere (cf. [5] 
or eq. above). This enabled us to define the EEP vortices as the points where 
the couplings are maximal, i.e., where the boundary dV^^^ of the quasi-Hermiticity 
domain intersects the circumscribed "upper-estimate" hypersurface dS^^\ In this 
context, a key technical result of paper I consisted in the derivation of the closed 
formula for all of the EEP coordinates (viz., of eq. ([H]) above). 

In the present continuation of paper I we intend to pay attention to the shape 
of the hypersurfaces dV^^^ in the vicinity of their EEP extremes. This is a well 
motivated project since, in spite of the reality (i.e., "mathematical" observability) 
of the energies of the system in its solvable EEP limit, the "physical" version of the 
same (i.e., strong-coupling) observability concept requires more than that. Obviously, 
during any measurement we should stay in the interior oi the domain V^^\ requiring 
that a small random perturbation of the couplings cannot induce a spontaneous 
complexification of some energies and a subsequent sudden collapse of the system. 

2.2 Observations made in paper |[6| 

Later on, we complemented the algebraic constructions of ref. [5] by the purely 
numerical study [6J of the possible complexification patterns of the spectra of our 
chain models We can summarize that 

• the perturbed harmonic-oscillator spectrum {En} remains all real in the weakly 
anharmonic regime characterized by a "sufficient" smallness of all the elements 
or couplings -ff^m+i at all m; 

• the pairs {En-, En+i) of energies coincide and complexify whenever their mutual 
coupling exceeds certain n— dependent "exceptional-point" value [T]. 

A particularly amusing empirical (though easily explained) observation made in [6] 
was that one can preserve the reality of the spectrum even when crossing the no- 
interaction boundary in accordance with the rule 

-yX, (A>0) (A<0). (13) 
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Although such a switch between the real and purely imaginary couplings will not be 
considered in what follows, we shall still use the parametrization fll3p in its strong- 
coupling real form iQ where A = 1 — t is assumed large while the new, formal 
auxiliary parameter t G (0, 1) remains, preferably, sufficiently small, characteriz- 
ing our PT— symmetric and up/down symmetric chain models H^'^^ in their most 
interesting, strongly non-Hermitian dynamical regime. 

In [6], the use of the parameters A or t has been shown to facilitate the study of 
the spectra, real or complex, at any A^. In addition, the use of the "renormalized" 
coupling strengths G„ opened the way towards a systematic (viz., combinatorial) 
classification of the non-equivalent energy-complexification patterns or, if you wish, 
of the non-equivalent scenarios of a "quantum catastrophe". Moreover, we found 
that once we fix a real J— plet of optional parameters (?„ in ansatz iQ, the range of 
the remaining free parameter t G (— oo, oo) splits in the four specific subintervals: 

1. in an "unobservable" regime, some of the energies are not real (so that H'^^^ 
itself is not quasi-Hermitian, QH) at t G oo, t{^QH){Gi, G2, . . . , Gj)^ ; 

2. the genuine quasi-Hermitian and PT— symmetric regime is encountered in the 
range of t G (t(Q/i-)(Gi, G2, • • • , Gj), t(PH)(Gi, G2, • • • , Gj)) where one stays 
safely inside T)^^\ The value of the parity-pseudo-Hermiticity boundary t(PH) 
is given by the decision [5] that the matrix H^^^ remains real, max„^„(t) < 1; 

3. in the next interval of t G (t(^PH){Gi,G2, ■ ■ ■ ,Gj),t(^H){Gi,G2, ■ ■ ■ ,Gj)^, the 
matrix ceases to be real and its V— (i.e., parity-) pseudo-Hermiticity 
"strengthens" to an 17— pseudo-Hermiticity. An A^ = 4 illustrative example of 
the "modified parity" r] (which may further vary with t) was given in ref. [6]; 

4. Hermitian regime enters the scene at t G (t(^H){Gi, G2, . . . , Gj), 00^. One has 
^n(t) > 1 at all n so that all the couplings become purely imaginary. 

In this setting, one notes a certain complementarity between the results of refs. [5] 
and [6]. In the former (and older) text we intended to stay safely inside the closure 
of the domain T)^^\ In this framework we showed that the choice of the optimal 
proportionality coefficients ([H]) in eq. f[T^ implies that we can minimize ^(qh) = 0. 
In this sense we found the maximal interaction strengths which are still compatible 
with the reality of the spectrum. In contrast, the scope of paper [6] was broader 
and covered (i.e., sampled, numerically) all the four eligible intervals of the auxiliary 
parameter t. 
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In what follows we shall show that and how the successful localization of the 
exceptional boundary point with t(Qiy) = in [5j can be extended to a certain 
approximative closed-formula description of all the ( J— 1)— dimensional boundary set 
in a vicinity of this point. This will definitely clarify, inter alii, that and how 
the necessary physical stability of our model can be guaranteed via a constructive, 
leading-order specification of an interior, strong- coupling part of the domain T)^^\ 

3 The method of exphcit constructions 
3.1 Guiding example: = 4 

In the first nontrivial N = 4 example we set gi = (1 — g2 = (1 — a) and 
consider 

^ -3 y/3-3P \ 

^(4)^ -1 2^1^ 

-2^/T^ 1 V3-3/^/ 

^ -V3-3/5 3 y 

with a , /? G (0, 1). This leads to the secular equation 

- (6^ + 4a) s- 36/? + 36a + 9/3^ = (14) 

with the doublet of available elementary roots, 

s± = 3P + 2a±2^3pa + a'^ + 9p-9a. (15) 

As long as the energies are square roots of these roots we must guarantee that s± > 0. 
Vice versa, the latter two inequalities may be understood as an implicit definition of 
the domain p(^) = V^*\a,p). More details may be found elsewhere [T5] . 

In an attempt to clarify the origin of ansatz (Q let us now introduce an auxiliary, 
redundant parameter t and set a = t a and (3 = tb, treating t as a radius of the EEP 
vicinity and preserving just the leading-order terms in t. Then, conditions s± > 
degenerate to the two elementary rules 

b>a + 0{t), a>b + 0{t) 

which may be interpreted as a requirement of a fine-tuned balance between a and b 
(or a and /?) near the EEP extreme. We may conclude that our ansatz is optimal 
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and that inside V^^'^ , the value of a can only differ from (3 in the next order in the 
small 

(3 = t + Bt^, a = t + At^. (16) 

Tractable as an exact change of variables (a,/?) — > {A,B), such a rule strictly 
replaces the inequalities s± > by the t— parametrized pair of the necessary and 
sufficient conditions of the physical acceptability of H^'^\ 

t^A^ + 5tA + 3t^AB + 4: + 3tB + 9B - 9A>0, (17) 

t^ + 2t^B + t^B^ -it'^B + 4:t'^A>0. (18) 

In the vicinity of the EEP extreme we can omit the higher-order terms from eqs. (fT7|) 
and (flSi) . This reduces this pair of inequalities to the compact leading-order estimate 

I>A-B>-1. (19) 

This leading-order rule defines V^^^ reliably near the EEP vertex. Indeed, in a graph- 
drawing scenario we may suppress the redundancy and fix t = /? (i.e., set B = 0). 
This converts f|T9l) into the explicit leading-order formulae which define the two 
branches of the boundary d'D^^\ 

9 4 

These two curves osculate at EEP so that the domain V^^^ appears sharply spiked 
near this extreme. 

It makes sense to believe that such a geometric property of the domain V is 
generic. At the higher matrix dimensions A^, precisely this hypothesis has been 
formalized by the tentative ansatz (Q. For a graphical illustration of its consequences 
we decided to sample the = 4 spectrum in Figure 1. At i? = it compares the 
quadruplet of the energies which are real at all t > (and which correspond to 
the "admissible" choice of the rescaled coupling A = 4/9 — 20/100 which lies safely 
inside T)^^^) with another quadruplet obtained at a "forbidden" A = 4/9 + 2/100 
[i.e., slightly outside violating the upper bound in eq. flTI?]) ] which remains all 
complex at all the positive small t < t(QH){A). 
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Figure 1: The dependence of energies, with two choices oi A ai N — A. 
3.2 Inessential changes at N = b 

At the first nontrivial odd dimension N — 5 with — b — 2 ^/l — /3, /3 e (0, 1) and 



g2 — a = JQ{1 — a), a e (0, 1), our model 



( -4 2v/r^ 

-2v/r^ -2 V6-6q; 

-V6-6q; V6-6q; 









V 













V6-6q; 2 2v^r^ 

-2^/r^ 4 



leads to the secular polynomial of the fifth degree in E which is divisible by E. Thus, 
one of the roots, viz., the energy i?2 = may be treated as trivial. The remaining four 
energies E may be computed from the two roots s — E"^ oi the quadratic equation 



(5), 



(21) 



where one easily evaluates 
^i^'H^i,^2) = 8/3 + 12q;, 



P2^{9\. = 48 q;/3 - 144 /3 + 144 a + 16 0^ 



We may skip the details ~ near EEP the construction of T)^^^ is strikingly analogous 
to its = 4 predecessor. 
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3.3 Inessential changes at = 6 

At = 6, three parameters a,/?,7 G (0, 1) enter the three couphng constants 

(71 = c = v^S (1 - 7) , g2 = h = 2p{l-(3), ^3 = a = 3v/l 
which specify the dynamics via the Hamiltonian 



a 
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(22) 



The shape of the associated domain of quasi-Hermiticity V = V{a, b, c) can be de- 
duced from the secular equation 



(23) 



det - E I^'^) =s'~3 Pf ) s' + 3 Pi') s - = , 
(notice the shghtly modified notation) re-written as the relation 

[s - si{a, b, c)] [s - S2{a, b, c)] [s - ^3(0, 6, c)] = 

between its roots Sk{a, b, c). As long as they define the energies E±k: = ±-y/ Sk{a, b, c), 
all of them must be nonnegative in T)^^\ For this reason, we have to satisfy three 
requirements Pj:^^ > 0, k = 1,2,3 plus a certain slightly more complicated fourth 
condition (with the derivation left to the reader as an exercise). 

From a geometric point of view, the construction of the boundary dV^^^ of the 
physical domain remains similar to its A^ = 4 predecessor. Thus, equation 

P^(^) = _ (^^2 ^ 2 6^ + 2 - 35) /3 = 

determines the ellipsoid which circumscribes the domain T)^^\ In contrast, the 
geometric interpretation of the further circumscribed surfaces (given by equations 
pi'"* = etc) is much less straightforward. For this reason we intend to employ 
ansatz IQ and to replace the variables a, j3, 7 by the new triplet A, B, C defined by 
the J = 3 version of the recipe, 

a = t + t^ + At^, f3 = t + t^ + Bt\ -i = t + t^ + Ct^. (24) 

The auxiliary variable t is redundant but useful because we intend to keep it small. 
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4 The method of rescaled couphngs 



Our above exact result (fT6|) can be reinterpreted as a tentative ansatz with the 
obvious generahzation (Q. The use of the latter rule requires an assumption of the 
smallness of |t| -C 1. This could simplify our secular eq. (ITOl) . near the EEP extremes 
at least. 

4.1 Guiding example: = 4 

In the new notation our secular equation (i.e., our implicit definition of the = 4 
energies) degenerates to the leading-order relation 

- 10 1 s + (36 A - 36 5 + 9) + 0{t^) = . 

After we introduce a new coupling parameter to = co^'^^ = 36 (A—B), this equation for 
the unknown quantity L = s/t acquires a transparent t— independent leading-order 
form which may be partially factorized, 

L2-10L + 9 + w = (L-l)(L-9) + cj = 0. (25) 

It is important that the roots of eq. fl25|) are known exactly, 

L± = 5 ± V16-UJ . 

4.1.1 The leading-order localization of the boundary dV^"^^ 

Obviously, the growth of > beyond its "upper limit" = 16 makes all the 
four roots L (i.e., all the related leading-order energies) complex. In the alternative 
scenario, the decrease of u < below its "lower limit" ujll = —9 makes just one of 
the roots (viz., L_) negative. 

Both these estimates coincide with the above-derived formula (fT9|) . We may 
conclude that the use of our perturbation-type ansatz (IT^ reproduces, completely, 
the leading-order information about the spiked shape of the boundary dV^^^ near the 
EEP extreme. 

We could also speak about the quadruplet of the t— dependent energies EQ(t), . . ., 
E^lt) studied in the strong-coupling regime. This dynamical regime is characterized 
by the small auxiliary quantities t which, in effect, measure the "distance" of our 
model H^^^ from its EEP t = extreme. In such an alternative language, the pair 



12 



of Ei{t) and E2(t) will complexify somewhere near u ~ uj^i, etc (cf. [6] for a more 
complete discussion and for a combinatorial classification of all the possible scenarios 
of complexification at any A^). 



4.1.2 A linearized localization of dV^^^ at small 



In the symmetric-coupling regime with A = B one stays safely inside the "physical" 
interior of T)^^\ Once we assume that the difference to ~ A—B is small, the extraction 
of all the energy roots becomes facilitated. In the leading-order approximation their 
a;— dependence remains linear, 

L_ = l + ^ + 0(a;2), L+ = 9-| + 0(c.2). (26) 

These formulae still offer a qualitatively correct perturbative explanation of the com- 
plexification pattern of the energies. Indeed, even from the oversimplified approxi- 
mation s = E'^ = Lt + 0{t^) using closed formulae (126!) one can deduce the rough 
estimates of the critical ujj^\ = 32 (or cca 23.4 in the second-order approximation 
in uj) and tu^^^ — ~^ (o^ —9.37 in the second-order approximation), yielding even a 
reasonably good quantitative prediction. 

4.1.3 A systematic build-up of higher-order corrections 

Whenever we keep the t— dependent version of our = 4 secular eq. (HM in full 
precision, we can construct the perturbation series for the energies at the small t in 
the standard manner, with 

oo oo 
k=0 p=0 

and with 



eP = V2 V95 -9A + 4 + 5 , 

Ei^^^^ ^ ' ( + +3B + 2A 

2^2 V95-9A + 4 + 5 W9B-9A + A j 

etc., or with 

Ei°^ = ^5-2^95"- 9 A + 4 , 



2^5-2V9B~-9A + i \ V9B-9A + 4J 

etc. These formulae offer another source of insight in the mechanisms of complexifi- 
cation of the spectrum at the boundary d'D^^\ 
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4.2 A change of the pattern at = 5 

The insertion of ansatz (1161) in secular equation (1211) gives 

) (6, a) = 20t + 8t^B + 12 t'^A . 

This confirms our expectations that the necessary non-negativity of this expression is 
guaranteed at all the not too large |y4|, \B\ ^ 1/t. Similarly we are able to evaluate 
and demand that, in full precision, 

{h, a) = 144 - 144 + 64 + 16 t'^B^ + 48 t^A + 48 t^A B + 80t^B >0. 

Obviously, the higher-order corrections may only be needed quite far from the EEP 
extreme. In the dominant order in t we may conclude that 

-20ts+ (144 A - 144 5 + 64) + 0{t^) = . 

Once we put s = Lt + (9(t^) again, we get 

L2-20L + 64 + e = (L-4)(L-16)+e = 0, e = lU {A - B) . (27) 

In the regime with a small e (of any sign), we easily get the leading-order e— dependence 
of both the roots. 

Marginally, it is amusing to notice that in the trivial case with A — B = e = the 
positive doublet of energies (as well as its negative counterpart) has a tendency of 
moving directly to the corresponding weak-coupling harmonic-oscillator limit. Our 
leading-order approximation E± = +^ (10 ± 6) t + 0(t'^) behaves as if being, para- 
doxically, exact at t = 1. 

4.2.1 Towards the boundary dV^^^ 

In a way paralleling our above = 4 considerations let us switch to a larger, 
positive e ^ A — B > 0. The two dominant energies > E2 > will then 
decrease and increase with the growth of e, respectively. One can predict an ultimate 
coincidence and a subsequent complexification of this doublet at a non-vanishing 
Kato's exceptional point localized at some "boundary coordinate" t = t(^QH) > of 
quasi-Hermiticity. 
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For the opposite, negative and decreasing £ < the approximative results are 
similar and do not lead to any contradictions. We may summarize that near the 
shared EEP maximum of gi and g2, the growth of the difference between a and /3 
would lead us out of the physical domain T>^^\ In the other words, the boundary 
dT>^^^ near EEP is of a sharply spiked form as well. 

5 Extrapolation hypothesis 

Let us now turn attention to the models with J = 3. We intend to offer some quan- 
titative arguments in favor of the intuitive idea that the exceptional-point bound- 
aries dV^^^ and dV^'^^ have a form of a surface of a deformed cube with protruded, 
razor-sharp edges (= double exceptional points) and with the spiked, strong-coupling 
vertices (i.e., triple exceptional points). 

We expect that such an intuitively transparent geometric interpretation of the 
shape of V^^^ and V^'^^ will provide a firm ground for extrapolations towards higher 
dimensions in subsequent sections of this paper. 

5.1 Expectable observations at N = 6 

After transition from = 4 to = 6, symbolic manipulations on a computer become 
rather lengthy. Still, they enable us to evaluate, in closed form, the approximate 
secular-equation coefficients 

) = )(c, b,a)^35t + 35 + Oit^) , 

P^^^ = )(c, b, a) = 259^2 + (216 ^ + 518 + 144 5 - 360 C) ^ + 0{t^) 

and 

Pg^^) = P^^\c, b, a) = (3600 A + 1200 C - 4800 B + 225) 

+ (3600 B - 1800 C + 675 - 1800 A) + 0{t^) . 
The resulting shortened leading-order secular equation is finally obtained and reads 

5^-35^5^ + 259^^5- (225 - oj^^^) = 0, oj^^'^ = 1200 (-3^ + 4P - C) . 

After we put s — E'^ — Lt + 0{t^), its partial factorization remains feasible, 

(L-l)(L-9)(L-25)+a; = 0. (28) 
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In a way paralleling the previous = 4 example, a transparency of this equation 
facilitates the analysis of the relationship between the variations of the coupling- 
dependent quantity uj ~ 45 — C — and of the physical spectrum of energies 
Eq = —E^ < El = —E4 < E2 = —E3, inside the domain of its reality at least. 

5.1.1 Towards the boundary dV^^^ 

The discussion is easier at the small to^^^ where all the three roots of eq. fl25]) are 
almost linear in u, 

8-24 ^ ^' 8-16 ^ ^' 16-24 ^ ^ 

At ~ and at the smallest t 7^ we have Eq ~ —5 ^/t, Ei ~ —3 y^, E2 ~ —Vt, 
E3 ~ Vi, £'4 ~ 3 and E^ 5 \/t. With the decrease of < the levels 
El and E2 (and also E^ and E^ get closer to each other and, in a way which 
parallels the similar observation made at = 4, they finally complexify at the 
critical uj^ll) ~ -323.1387184 near E^ ^ E^ ^ 2.147400716 v^. In the light of 
the obvious fact that the growth of every coupling causes an attraction of the 
corresponding levels [9], the interpretation of the above rule is easy and consistent 
because the decrease of uj means not only a smaller B < (3 A + C)/4 but also, at 
any given t, gi and g^, an enhancement of the coupling g2 between the levels which 
complexified. 

In the opposite direction, the increase of > causes the decrease of the small- 
est root Ls{uj) (= El^/t) of eq. fl28l) while the other two roots Lm{uj) (= Ef^/t) 
and Ld{uj) (= Eq^/I) start approaching each other. In this way, the first and 
dominating complexification involves the pair £^2,3 (with £^2,3 ~ 0) because it takes 
place exactly at uJ{ul) = 225, i.e., much earlier than the subsequent ultimate com- 
pletion of all the complexification process with E^ ^ Ei ^ 4.326893054 near 
t^{c/L[4,5]) ~ 1081.657240. 

5.2 Expectable changes at A^^ = 7 

With N = 7 and with the J = 3 rule ([MD, we employ the symbolic manipulations 
on a computer and get five terms in 

pf ) (a, 6, c) = 56 1 + 56 + 0{t^) 
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sixteen terms in 

pP{a,b,c) = 784^2 + 

and 33 terms in 

pP (a, b, c) = (7200 C - 21600 B + 14400 A + 2304) + 0{t^) . 
The resulting shortened leading-order secular equation reads 

- 56 1 + 784 s - (2304 -e)t^ = 0, e = 7200 {-C + 3 B - 2 A) . 
With s = L t it factorizes as follows, 

(L-4)(L-16)(L-36)+£ = 0. (29) 

5.2.1 Towards the boundary dV^'^^ 

The three closed solutions of eq. (129]) may be expanded in the powers of e, 

384 ^ ^' 240 ^ ^' 640 ^ ^ 

A balanced return to the unperturbed values occurs now along a "middle line" with 
e = 0, i.e., for B = {2A + C)/3. 

For the diminished -Bs (i.e., for a stronger coupling b = g2 between Ei and E2), 
we are getting closer to the boundary surface along which the levels Ei and E2 merge 
and subsequently complexify. This behaviour is confirmed by our formula because 
the shift e becomes negative in such a scenario. 

The growth of e > 0, on the contrary, may be assigned to the diminished A or C. 
In the former case the growth oi a = gi implies that Ei —* while the alternative 
growth of c = results, naturally and expectedly, in the merger of E2 with E3. 

6 Arbitrary even dimension and the localization 
of boundaries dV^^"^^ in the EEP vicinity 

We saw that in many a respect one should pay separate attention to the models with 
even and odd A^. Thus, we shall search for a successful extrapolation of our low— 
"experiments" in two separate sections, starting with the family where A^ = 2 J. 
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Our first assumption concerns the values of tlie rescaled couplings (?„ (abbrevi- 
ated, occasionally, as A, B etc) which will be assumed bounded while the scahng 
parameter t itself will be assumed small. 

Our second, main assumption is that our older eqs. fl25l) and eq. fl28|) are just the 
respective iV = 4 and at = 6 special cases of the general leading-order secular 
equation 

(L- 52) ... (l-[2J-1]2) + cu(-^) = 0. (30) 

We tested this hypothesis together with the ansatz ^ by extensive symbolic manipu- 
lations which confirmed its validity at a// the integers = 2 J in a way documented in 
Table 1 . This Table summarizes also the corresponding resulting explicit formulae for 
the J— dependent quantities uj'^'-'\ Using the same extrapolation philosophy as above, 
we may extract the universal algebraic definition of the coupling-characterizing ex- 
pression uj^'^^ as given by the following extrapolated conjecture at all the positive 
integers J, 

u^'^ = 2 (2J - 1) (2J - 1)! f (-1)^-+! [ 2 - 2 \ ^^^^ ^3^^ 

n=i \J -2 + n I 

This is our first main result deduced via an interpolation of the computed coefficients 
at the first few dimensions A^, complemented by a subsequent (and, of course, much 
easier) verification of the hypothesis at a number of the higher A^. Formula (!3T|) 
contains just the usual combinatorial numbers (given by the standard Pascal-triangle 
recurrences) and the anomalous scaling factor 9{n) which is equal to one at ri > 1 
and to one half at n = 1. 

The geometric interpretation of formula fl3T]) could be discussed in an immediate 
parallel with the texts on the special cases with A^ = 4 and A^ = 6 in the respective 
paragraphs 14.11 and 15.11 At the general J, a careful inspection of the first few sec- 
ular eqs. fISUl) seems to indicate that at the smallest ts, the complexification of the 
spectrum always proceeds 

• through the single merger, at the zero energy, of the "middle" levels 
and Ej (at some > for J = 1, 3, 5, . . . and at some ^fii) < ^ 
J = 2,4,6,...), 

• through the two simultaneous, symmetric mergers of the negative £"7-3 with 
Ej^2 and of the positive -Ej+i with -E'j+2 (at < for J = 1, 3, 5, . . . and 
at uj\1j\. > for J = 2,4,6,. . .). 



18 



7 All odd dimensions and the boundaries 
in the EEP vicinity 

When the dimension is odd, = 2J + 1, the apphcation of the ansatz (Q sim- 
phfies the problem significantly as well. This ansatz can be perceived again as an 
equivalence transformation based on a mere change of the variables, Qn Gn- It 
maps the domain T)^^\gi, . . . ,gj) of quasi-Hermiticity oi H = H^'^'^^^\gi^ . . . ,gj) 
into another, equivalent manifold V^'^'^'^^^Gi, . . . , Gj) of the acceptable parameters 
mH = H^''+'\G^,...,Gj). 

For all the odd dimensions we may ignore the persistent "middle" bound-state 
energy Ef"^^^^ = 0. Mutatis mutandis, this enables us to use the same J and to apply 
the same (or at least very similar) geometric, algebraic and analytic considerations 
as above. 

One of the most visible differences between the models if^^-^) and may 
be traced to the fact that in place of the (hyper)ellipsoids pertaining to the even 
dimensions we have to deal with the simpler (hyper) spheres S^^^ at all = 2J + 1. 
In the same comparison, in contrast, the leading-order secular equation 

(l - 2^) (l - 42) (l - 6^) . . . (l - [2 J]2) + = (32) 

pertaining to the odd = 2J + 1 is slightly more complicated. This is demonstrated 
by Table 2 showing that the general formula for the factors e*-'^-* 

eW = 2(2J-l)(2J)! X:(-1)"""C^;5G'„ (33) 

n=l 

is perceivably less explicit. The reason is that the elementary combinatorial coeffi- 
cients in eq. ( 13T|) are replaced now by no-name integer coefficients C(^)- 

Fortunately, the latter coefficients can be understood as an immediate general- 
ization of the current combinatorial numbers since, firstly, they exhibit a left-right 
antisymmetry C*(/Z„) = ~G^^j at = 1, 2, . . . and, secondly, these coefficients may 
be generated by the recurrences 

from the initial row cj^^j at J = 1 which contains a unit at = 1 and zero for n > 1 
(cf. Table 3). Thirdly, a decomposition of the sum flMl) in the two steps specifies 
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the coefficients C^^^ as a subset of the triangle generated by the usual combinatorial 
three term recurrences. The latter generation pattern is sampled in Table 4. One 
notices that both the three-term recurrences and their illustrative Table 4 differ 
from their standard Pascal-triangle predecessors merely in an anomalous, left-right 
antisymmetric initialization. 

Once we return to our leading-order energies E^'^^^^ we may say, on the basis of 
eq. ([22D, that their /eadm^'- orc?er complexification always seems to proceed in one of 
the following two ways: 

• through the two simultaneous, symmetric mergers of the negative -Ej_2 with 
Ej-i and of the positive -Ej+i with -E'j+2 (at < for J = 1, 3, 5, . . . and 
at ej;^^^) > for J = 2,4,6,...). 

• through the two simultaneous, symmetric mergers of the negative -Ej_3 with 
Ej_2 and of the positive -Ej+2 with -Ej+3 (at > for J = 1, 3, 5, . . . and 
at ^l^l-) < for J = 2,4,6,...). 

In the other words, the implicit definition (!32|) of the spectrum has a similar geometric 
interpretation and implications for the shape of the boundary dT>^'^'^^^'> as its even- 
dimensional predecessor (!30!) did. 

8 Discussion 

8.1 The real Kato's exceptional points 

In the Kato's book [1], certain k— dependent one-parametric families of by 
matrices were discussed, with elements defined as holomorphic functions of k in a 
complex domain Dq. It has been noticed there that up to not too many "exceptional 
points" k'^^^'> E Dq, the number of eigenvalues of any such matrix H{k) is equal to 
a constant. 

One of the Kato's illustrative two-by-two examples possessing just a real pair 
of the Kato's exceptional points was discussed in our recent paper [H]. In his/our 
example (= H^'^'^ in our present notation) we worked just with the real k = gi and 
identified the Kato's exceptional points as points of the boundary of the correspond- 
ing quasi- Hermiticity domain, k^^^^ G dV^'^\ 
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This was one of the important motivations of our continuing interest in the struc- 
ture of the sets dV^^^ of the real exceptional points attached, in particular, to the 
real chain models H^^^ at the higher matrix dimensions A^. In these particular mod- 
els, the sets of the real Kato's exceptional points (i.e., the boundary sets dV of the 
models in question) exhibit a nice and transparent hierarchical pattern of confluence. 

Among all the multiple exceptional points, an extreme is represented by the 
J— tuple confluence of these points at the strong-coupling extremal EEP vortices 
[5]. Of course, all the similar multiple confiuence(s) of the exceptional-point hy- 
persurfaces of various dimensions are in a close correspondence with some physical 
critical phenomena. We might emphasize that in our family H^'^\ a fine-tuning 
mechanism emerges which aligns the values of the physical couplings in a manifestly 
non-perturbative though still analytically tractable and presumably also generic and, 
in its qualitative aspects, not too model-dependent manner. 

8.2 Separable anharmonicities 

On the more practical side of our present considerations and constructions it is 
worth emphasizing that the Kato's abstract models of a generic dependence of the 
eigenvalues on a single variable parameter in if = H{k) ^ W need not always offer 
a better qualitative understanding of the situation. 

For an explicit illustrative example we may return to the differential-equation 
example of paper [2]. An inspection of the Table 1 of this paper reveals that once 
we vary just an arbitrarily chosen free parameter a (as defined in eq. ([3]) above), no 
obvious pattern is detected in the complexifications and decomplexifications of the 
energy levels which occur in the Table in an apparently unpredictable and more or 
less chaotic manner. 

What we would recommend in similar situations would be, in the light of our 
present experience, a deeper inspection of the separate matrix elements themselves, 
followed by a subsequent tentative re-definition of the spectrum as depending on some 
more (say, J) independently variable parameters m. H = H{ki, K2, . . . , nj) 7^ H^. 
Only then one would have a right to expect that one finds an "optimal" linear 
combination of kjS which wouldn't cross the global boundary &D at random or, at 
least, which would stay locally more or less perpendiculat to dV. 

One could expect the existence of an "optimal" size J of the set of the auxil- 
iary parameters in many other manifestly non-Hermitian models based on the use 
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of a complex potential V{x). Unfortunately, it is difficult to predict how such a 
type of analysis would prove efficient in practice. In this sense a return could be 
recommended to many older papers on the subject [16]. 

The problem looks interesting even as a purely mathematical puzzle of survival of 
the reality of the energies (i.e., in principle, of the observability of the system) after 
real potentials are replaced by their complex PT— symmetric versions. An interest 
in such an apparent paradox has been evoked almost ten years ago [121 [H] but it 
took several years before the necessary rigorous proofs of the reality of the spectrum 
became available [T8] . 

Bad news is that too often, the reality proofs are fairly complicated [19j. In 
parallel, good news is that once all the energies stay real, one returns easily to the 
standard principles and formalism of Quantum Mechanics by using an ad hoc scalar 
product in the physical Hilbert space 7i. This means that one merely replaces the 
usual overlaps (0 | tp) by their generalizations {(f)\Q\ilj) (this idea belongs to Scholtz 
et al [n]). 

Although, as a rule, the physical Hilbert-space metric = 0^ > proves compli- 
cated and non-local, people often succeed in its (e.g., perturbative [20j) construction. 
Of course, the explicit specification of the whole "physical" domain V seems to be 
an even more difficult task. That's why we recommended here the use of the finite- 
dimensional prototype matrices H^^\ 

It is obvious that our knowledge of the boundaries dT> = dT>{H) of the physical 
consistence of a generic "PT— symmetric Hamiltonian is necessary for any reliable 
physical prediction or for a consistent probabilistic interpretation of the results of 
quantum measurements. This boundary may have a fairly complicated shape even 
in the simplest, exactly solvable models [21]. At the same time we believe that 
the present constructive clarification of the geometric structure of the corresponding 
"prototype physical horizons" dV^^^ {jl^^^^ (where even the diagonalizability of the 
matrices H^^'^ themselves gets lost!) could accelerate the progress in the analysis of 
the more common Hs defined via potentials. 

9 Summary 

Without any limitations imposed upon A^, the series of all the models H^^^ has 
been shown to admit an innovative strong-coupling approximative treatment. In 



22 



particular, the underlying polynomial secular equations were shown to degenerate to 
the closed leading-order form in the strong-coupling dynamical regime. 

An appropriate technical tool for our simplification of the chain-model bound- 
state problem in question has been found in ad hoc perturbation ansatz iQ which 
reparametrizes all the coupling constants in a way inspired by their available closed- 
form t —>■ limit known from our previous paper [H]. 

Also our present methods were inspired by the same reference - we further 
developed the application of the brute-force symbolic-manipulation techniques as 
well as of the systematic formulations and verifications of extrapolation hypothe- 
ses. Marginally, we may add that one of our results, viz., the implicit leading-order 
formula for the energies in odd dimensions N = 2J + 1 required a slightly unusual 
though straightforward and immediate generalization of the standard binomial coef- 
ficients. Our new combinatorial coefficients were shown closely related to a pair of 
alternative generalizations of the popular Pascal triangle. 

Via a deeper analysis of the leading-order secular polynomials we were, subse- 
quently, able to demonstrate, in the whole strong-coupling dynamical regime 
that and how the determination of the asymptotic strong-coupling parts of the do- 
mains of quasi- Hermit icity T)^^^ degenerates to elementary approximative formulae 
at all the dimensions A^. 

Our main conclusion is that our class of models H^^'^ could be understood, in 
many a respect, as representing or mimicking many generic features of the general 
PT— symmetric and pseudo-Hermitian Hamiltonians. In this sense, in particular, our 
present quantitative description of the mechanism of complexification of the energies 
could find an important application as a guide to our understanding of the change 
of the dynamical regime called the spontaneous breakdown of PT— symmetry [21] . 
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Table captions 

Table 1. Auxiliary functions of couplings cj*^^^ (even N = 2J). 
Table 2. Auxiliary functions of couplings e^^^ (odd N = 2J+1). 
Table 3. Coefficients C^f^^ of eqs. ([MD and (l34l). 
Table 4. Pascal-like triangle, with coefficients C^^^j* underlined 

Figure captions 

Figure 1. The t— dependence of energies, with two choices of 
^ at iV = 4. 



24 



References 

[1] T. Kato, Perturbation Theory for linear Operators (Springer, Berlin, 1966). 

[2] Y. Ben-Aryeh and R. Barak, Phys. Lett. A 351 (2006) 388. 

[3] C. M. Bender, Reports on Progress in Physics 70 (2007) 947; 

J. Phys. A: Math. Gen. 39 (2006), Nr. 32 (ed. H. Geyer, D. Heiss and M. Znojil); 
Czechosl. J. Phys. 56 (2006), Nr. 9 (ed. M. Znojil). 

[4] A. Mostafazadeh, J. Math. Phys. 43 (2002) 205 and 2814 and 3944; 
L. Solombrino, J. Math. Phys. 43 (2002) 5439. 

[5] M. Znojil, J. Phys. A: Math. Theor. 40 (2007) 4863. 

[6] M. Znojil, Phys. Lett. B 650 (2007) 440. 

[7] J. H. Wilkinson, The algebraic eigenvalue problem (Clarendon, Oxford, 1965). 

[8] C. M. Bender, D. C. Brody and H. F. Jones, Phys. Rev. Lett. 89 (2002) 0270401; 
Q.-H. Wang, Czechosl. J. Phys. 54 (2004) 143; 
A. Mostafazadeh and S. Ozcelik, Turk. J. Phys. 30 (2006) 437. 

[9] M. Znojil and H. B. Geyer, Phys. Lett. B 640 (2006) 52 and 

M. Znojil and H. B. Geyer, Phys. Lett. B 649 (2007) 494 (erratum). 

[10] M. Znojil, Phys. Lett. B 647 (2007) 225. 

[11] F. G. Scholtz, H. B. Geyer and F. J. W. Hahne, Ann. Phys. (NY) 213 (1992) 
74. 

[12] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80 (1998) 5243. 

[13] O. Kirillov and U. Giinther, J. Phys. A: Math. Gen. 39 (2006) 10057. 

[14] D. Heiss, J. Phys. A: Math. Gen. 39 (2006) 10077. 

[15] M. Znojil, Phys. Lett. A 367 (2007) 300. 



25 



[16] E. Caliceti, S. Graffi and M. Maioli, Commun. Math. Phys. 75 (1980) 51; 
V. Buslaev and V. Grecchi, J. Phys. A: Math. Gen. 26 (1993) 5541; 
G. Alvarez, J. Phys. A: Math. Gen. 27 (1995) 4589; 

F. M. Fernandez, R. Guardiola, J. Ros and M. Znojil, J. Phys. A: Math. Gen. 
31 (1998) 10105; 

E. Delabaere and F. Pham, Phys. Letters A 250 (1998) 25. 

[17] C. M. Bender and K. A. Milton, Phys. Rev. D 55 (1997) R3255; 

C. M. Bender, S. Boettcher and P. N. Meisinger, J. Math. Phys. 40 (1999) 2201. 

[18] P. Dorey, C. Dunning and R. Tateo, J. Phys. A: Math. Gen. 34 (2001) L391. 

[19] P. Dorey, C. Dunning and R. Tateo, J. Phys. A: Math. Gen. 34 (2001) 5679; 
K. C. Shin, Commun. Math. Phys. 229 (2002) 543. 

[20] A. Mostafazadeh and A. Batal, J. Phys. A: Math. Gen. 37 (2004) 11645; 

A. Mostafazadeh, J. Phys. A: Math. Gen. 39 (2006) 13495; 

T. Curtright and A. Veitia, J. Math. Phys., to appear 
(|arXiv:q uant-ph/ 070T006|y2). 

[21] G. Levai and M. Znojil, Mod. Phys. Lett. A 16 (2001) 1973; 

M. Znojil and G. Levai Mod. Phys. Lett. A 16 (2001) 2273; 

A. Khare and U. Sukhatme, J. Math. Phys. 46 (2005) 082106; 

P. Dorey, A. Milhcan-Slater and R. Tateo, J. Phys. A: Math. Gen. 38 (2005) 
1305; 

A. Sinha and P. Roy, J. Phys. A: Math. Gen. 39 (2006) L377; 
A. Fring, Mod. Phys. Lett. A 21 (2006) 691. 



26 



Table 1: Auxiliary functions of couplings uj^'-'^ (even = 2 J). 



J 


dimension 


cu(-^)/[2(2J-l)(2J-l)!] 


2 


4 


A-B 


3 


6 


-3 A+4 B-C 


4 


8 


10A-15B+6C-D 


5 


10 


-35A+56B-28C+8D-E 


6 


12 


126 A-210 B+120 C-45 D+10 E-F 



Table 2: Auxiliary functions of couplings e'-"'^ (odd A^ = 2 J + 1). 



J 


dimension 


£(^)/[2(2J-l)(2J)!] 


2 


5 


A-B 


3 


7 


-2A+3B-C 


4 


9 


5A-9B+5C-D 


5 


11 


-14A+28B-20C+7D-E 


6 


13 


42 A-90 B+75 C-35 D+9 E-F 



Table 3: Coefficients Crl of eqs. (|33|) and (l3il). 





n= 


... -3 


-2 


-1 





1 


2 


3 


4 5 ... 


J 






(1) 






(0) 


(0) 


(-1) 


(1) 


(0) 


(0) 




2 








-1 


-1 


1 


1 






3 






-1 


-3 


-2 


2 


3 


1 




4 




-1 


-5 


-9 


-5 


5 


9 


5 


1 


5 




-1 -7 


-20 


-28 


-14 


14 


28 


20 


7 1 
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Table 4: Pascal-like triangle, with coefficients C}J underlined 



J 






n 




1 




2 


3 


(1) 




-1 


-1 





1 


1 






2 


-1 




-1 




1 




1 






-1 


-2 









2 


1 




3 


-1 -3 




-2 




2 




3 


1 




-1 -4 


-5 









5 


4 


1 


4 


-1 -5 -9 




-5 




5 




9 


5 1 
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